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Abstract 

An  apparatus  is  described  which  has  been  designed  and  built  to 
enable  measurement  of  the  thermal  conductance  of  soils  at  hot  face 
temperatures  approaching  1700  °C.  The  method  utilizes  radial  heat 
flow  through  a hollow  cylinder  of  soil  contained  between  a central 
ceramic  core  and  an  outer  water-cooled  metal  shell.  Apparent  thermal 
conductivity  values  are  reported  for  nine  natural  or  artificial  soils 
representative  of  subsoils  found  on  most  of  the  earth's  surface.  For 
four  of  these  soils,  measurements  have  been  completed  on  two  different 
specimens . 

Keywords:  Conductance,  conductivity,  heat  transfer,  nuclear 

safety,  soils,  temperature,  thermal  conductance, 
thermal  conductivity. 


1.  Introduction 


In  conjunction  with  safety  evaluation  of  space  nuclear  power  systems, 
it  is  desired  to  be  able  to  predict  the  maximum  steady-state  temperature 
which  will  be  reached  by  a given  nuclear  power  supply  after  reentry 
impact  burial  in  earth.  In  order  to  do  this  it  is  necessary  to  have 
information  regarding  the  thermal  conductance  of  soils  under  conditions 
in  which  heat  flows  from  a source  which  may  be  as  hot  as  1700  °C  to  a 
sink  at  ambient  temperatures, 

2.  Purpose  and  Scope 

Contract  No,  ASB48-6179  encompasses  an  experimental  study  to  measure 
the  effective  thermal  conductance,  as  a function  of  temperature,  of  nine 
soils  (natural  and  artificial)  selected  and  furnished  by  Sandia  Corpora- 
tion, Test  temperatures  shall  be  from  room  temperature  to  just  below  the 
molten  range  for  each  soil,  with  a maximum  as  near  1700  °C  as  possible 
wherever  the  melt  range  level  permits.  Duplicate  runs  shall  be  made  on 
each  type  of  soil  using  different  samples. 

The  thermal  conductance  shall  be  measured  in  a manner  that  simulates 
a heat  source  buried  in  soil.  Artificial  gas  pressure  or  inert  gas  shall 
not  be  used  in  these  measurements.  Application  of  any  compaction  load  to 
the  soil  samples  during  test  is  not  allowed. 


3,  Soil  Samples 


Samples  of  the  following  nine  soils  selected  by  Sandia  Corporation 
have  been  received  at  the  National  Bureau  of  Standards: 

1.  Calcareous  Soil  (natural  weathered  limestone) 

2.  Granitic  Detrital  Soil  (weathered  decomposed  granite  soil) 

3.  Dune  Sand  (windblown  sand) 

4.  Magnesian  Soil  (magnesium  aluminum  silicate) 

5.  Podzol  Soil  (leached  organic  timberland  soil) 

6.  Coastal  Plains  Clay  (ccastal  flood  plain  soil) 

7.  Laterite  Soil  (tropical  rain  forest  soil) 

8.  Estancia  Playa  (Dog  Lake)  Soil  (highly  saline  playa  soil) 

9.  Ottawa  Sand  (silica-artificial  soil) 

All  of  these  soils  were  screened  to  a maximum  particle  diameter  of  1/16 
inch  (0.16  cm)  prior  to  shipment  to  NBS. 

Sandia  has  provided  us  with  detailed  information  concerning  each 
soil  sample.  These  data  include  theoretical  estimates  of  solid-state 
reaction  temperatures  and  reaction  magnitudes  and  of  the  melt  range  for 
each  soil  to  be  tested.  In  order  to  experimentally  investigate  the  soil 
melt  ranges,  significant  reaction  temperatures,  and  rates  and  magnitudes 
of  reactions,  differential  thermal  analyses  (DTA)  of  each  soil  sample 
were  carried  out  by  the  Analytical  Coordination  Chemistry  Section  of  the 
NBS  Analytical  Chemistry  Division,  The  soil  sample  characterizations 
supplied  by  Sandia  Corporation  and  the  results  of  the  NBS  DTA  determina- 
tions are  given  in  Appendices  A and  B. 
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4.  Method  and  Apparatus  for  Thermal  Conductance  Measurements 

4,1.  Philosophy  of  Approach 

Our  prime  concern  in  deciding  on  the  best  experimental  approach  to 
this  problem  was  to  simulate  as  closely  as  possible  the  conditions  which 
a buried  heat  source  would  experience.  The  following  considerations 
apply: 

a.  The  maximum  steady-state  temperature,  (°C),  at  the  surface 
of  a heat  source  which  is  surrounded  by  soil  or  other  material 
is  given  by 


F£ 

^ - Ts  + 


(1) 


where  Tg  (°C)  is  the  temperature  of  the  sink  to  which  heat  is 
lost.  The  quantity  Q (W/m2)  is  the  power  output  per  unit  area 
from  the  source.  The  geometrical  factor,  F,  depends  on  the 
spatial  configuration  of  the  heat  source  and  sink.  The  mean 
thermal  conductivity.  A*  (W  m”*  °C“^) , is  defined  by 


K*  - \*  (T  , T ) 

8 m 


1 T 

— ~ T ^ MDdT, 


(2) 


ms  ■ 

where  A(T)  is  the  temperature-dependent  thermal  conductivity  of 
the  soil.  From  this  we  see  that  it  is  not  necessary  to  measure 
the  actual  thermal  conductivity,  \(T),  as  a function  of  tempera- 
ture if  the  mean  thermal  conductivity,  A*(Tg,Tj|),  can  be  measured 
directly  as  a function  of  the  hot  side  temperature.  This  can  be 
done  by  measuring  Tm  as  a function  of  Q for  some  configuration 
having  a known  geometrical  factor,  F,  and  then  computing  \ 
using  eq„  (1). 
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b.  If  X rat^er  than  X(T)»  is  measured,  the  cold  surface 

of  the  specimen  can  be  kept  at  room  temperature  thus  consider- 
ably simplifying  the  problem  of  containment. 

c0  If  X(T)  is  desired,  it  can  be  computed  from  measured  values  of 

**<Ts.V- 

d.  Since  the  quantity  of  prime  interest  is  the  maximum  temperature 

■u* 

attained  by  the  buried  radio-active  source,  the  values  of  \ 
which  are  appropriate  are  those  obtained  after  monotonic  heating 
of  the  hot  surface  of  the  specimen. 

e.  The  specimen  size  should  be  large  enough  that  the  heat  flow 
path  is  much  larger  than  the  average  particle  size  of  the  soil. 
For  a maximum  particle  diameter  of  1/16  inch  (0.16  cm),  it  would 
seem  appropriate  for  the  sample  to  be  at  least  1 inch  (2.5  cm) 
thick  from  the  hot  surface  to  the  cold  surface.  Similarly,  the 
heat  source  in  contact  with  the  specimen  should  be  large  in 
comparison  with  the  average  particle  size  of  the  soil.  This 
would  preclude  the  use  of  "hot-wire*'  methods. 

f.  Since  the  final  steady-state  temperature  at  the  source  surface 
is  of  principal  concern,  it  would  not  be  appropriate  to  attempt 
to  derive  thermal  conductivity  values  from  measured  values  of 
thermal  diffusivity  and  specific  heat.  This  is  particularly  so 
for  materials  such  as  soils  which  are  subject  to  decomposition, 
phase  changes,  moisture  migration,  and  other  phenomena  involving 
heat  generation  or  absorption. 
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g.  However,  it  must  be  realized  that  a buried  source  may  experience 
transient  temperatures  which  are  higher  than  the  final  steady- 
state  temperature  at  the  source  after  equilibrium  is  attained,, 

This  is  so  because  the  thermal  conductivity  may  increase,  e.g., 
due  to  sintering,  with  time  at  any  given  temperature.  The  best 
way  to  determine  what  would  be  the  maximum  transient  temperature 
rise  at  the  surface  of  a buried  heat  source  of  strength  per  unit 
area,  Q(W/m"  ) , would  be  by  direct  measurement  of  the  temperature- 
f time  relationship  in  an  experiment  which  closely  simulates  the 

conditions  which  a heat  source  buried  in  soil  would  encounter. 

4.2.  Method  & General  Mathematical  Analysis 
It  was  decided  to  measure  the  thermal  conductance  of  the  soil  samples 
using  the  method  of  radial  heat  flow  in  a right  circular  cylinder.  A 
cross  section  of  the  apparatus  is  shown  in  figure  1.  The  specimen  is  con- 
tained within  the  annular  space  between  the  outer  radius,  a,  of  a ceramic 
core  and  the  inner  radius,  b,  of  a brass  shell.  A measured  quantity  of 
heat  per  unit  time  generated  electrically  in  the  ceramic  core  flows  radially 
through  the  specimen  to  the  inner  concentric  water-cooled  brass  shell. 

The  ceramic  core  has  a concentric  ring  of  equally  spaced  holes  at  a radius, 
r',  parallel  to  the  axis,  each  containing  a heater  wire.  Temperatures 
are  measured  by  an  axial  thermocouple  in  the  ceramic  core  and  by  thermo- 
couples attached  to  the  outer  surface  of  the  inner  brass  shell  at  radius  c. 
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Figure  1.  Horizontal  cross  section  of  the . apparatus  for  measuring 
the  thermal  conductance  of  soils. 
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In  figure  1,  if  the  cylindrical  surfaces  r * a and  r = b are  each 
isothermal,  the  heat  flow  between  these  surfaces  will  be  radial  except 
near  the  ends  of  the  apparatus.  In  general,  the  thermal  conductivity 
of  the  specimen  material  will  vary  with  temperature,  and  the  heat  flow 
rate  per  unit  area  through  a cylindrical  element  of  the  sample  is 


where  temperature  is  denoted  by  the  symbol  T and  the  temperature-dependent 
thermal  conductivity  by  x(T).  If  eq.  (3)  is  rearranged  and  integrated, 
we  obtain 


where  and  are  the  temperatures  at  r * a and  r - b respectively. 


(3) 


X (T)dT  - QaMb/a)  , 


(4) 


Let  us  define  a mean  thermal  conductivity  \ (T|j,Ta)  over  the  temperature 
range  to  Ta : 


X (T)dT 


(5) 


Combining  eqs.  (4)  and  (5),  we  obtain 


(6) 


«n« 

as  the  expression  which  can  be  used  to  compute  \ from  measurable 


quantities . 
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As  an  example  of  the  use  of  (6)  let  us  consider  the  case  where  the 

thermal  conductivity  of  the  sample  can  be  assumed  to  vary  linearly  with 


temperature: 


X(T)  = XoCl  + Cfc(T  - T0)]  , 


(7) 


where  xG  *s  the  thermal  conductivity  at  an  arbitrary  reference  tempera- 
ture, T0,  and  0^  is  the  corresponding  temperature  coefficient  of  thermal 
conductivity.  Substitution  of  (7)  into  (5)  yields,  after  integration. 


Thus,  for  a material  for  which  the  thermal  conductivity  varies  linearly 
with  temperature,  eq„  (6)  yields  a thermal  conductivity  corresponding  to 
the  average  of  the  temperatures  at  the  hot  and  cold  surfaces  of  the 


Equation  (6)  requires  a knowledge  of  the  temperatures  and  Tb  at 
the  inner  and  outer  surfaces  of  the  sample.  Experimental  difficulties 
(e.g.,  finite  size  of  temperature  sensors  and  contamination  of  sensors 
by  the  sample)  preclude  direct  measurement  of  Tfl  and  Tb.  Reference  to 
figure  1 shows  that  there  are  four  temperature  drops  which  must  be  con- 
sidered in  deriving  Ta  and  from  the  measured  temperatures: 

1.  The  temperature  drop  between  the  thermocouple  well  in  the  center 
of  the  ceramic  core  and  the  surface,  r * a,  of  the  ceramic  core. 


X (Tb,Ta)  ” X 0[  1 + Qo(' 


Ta  + Tb 


2 


sample. 
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2„  The  temperature  drop  due  to  the  thermal  contact  resistance  between 
the  surface  of  the  ceramic  core  and  the  inner  surface  of  the  sam- 
ple, both  surfaces  being  nominally  at  r = a. 

3.  The  temperature  drop  due  to  the  thermal  contact  resistance 

between  the  outer  surface  of  the  sample  and  the  inner  surface 
of  the  brass  shell,  both  surfaces  being  nominally  at  r * b, 

4„  The  temperature  drop  between  the  inner  surface,  r = b,  of  the 

brass  shell,  and  the  outer  surface,  r = c,  where  the  tempera tur? 
is  measured. 

If  the  circle  of  heater  wires  at  r = r * were  a continuous  cylindrical 
heat  source,  the  entire  region  r < r*  inside  the  heater  circle  would  be 
isothermal.  For  the  case  of  a finite  number  of  line  heat  sources  at 
r = r 1 , the  temperature  measured  at  the  axis  is  equal  to  the  average 
temperature,  T^,  at  the  radius  of  the  heater  circle  [l,2]1.  Therefore  the 
temperature  drop  between  the  thermocouple  well  in  the  center  of  the  core 
and  the  surface  of  the  core  is  given  by 


Th  - V - k 


Oaln(a/r ') 
c 


> (9) 

where  kc  is  the  thermal  conductivity  of  the  ceramic  core  material  and  Ta 1 
is  the  temperature  at  the  outer  surface,  r = a,  of  the  ceramic  core  (as 
opposed  to  Ta,  which  is  the  temperature  at  the  inner  surface,  also  r 35  a, 
of  the  sample). 


1 Figures  in  brackets  indicate  the  literature  references  at  the  end 
of  this  report. 
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If  we  designate  R as  the  thermal  contact  resistance  per  unit  area 

3 

( € nf  at  the  core  specimen  interface  (r  = a),  then  the  temperature 

drop  across  this  interface  is  given  by 

Ta ' - Ta  " RaQ  . (10) 

Similarly,  the  temperature  drop  at  the  specimen* she 11  interface  (r  ■ b) 

is  given  by 

Tb  - V “ f RbQ  , (U) 


where  T^*  is  the  temperature  at  the  inner  surface,  r * b,  of  the  brass 
shell  (as  opposed  to  T^,  which  is  the  temperature  at  the  outer  surface, 
also  r = b,  of  the  sample)  and  is  the  thermal  contact  resistance  per 
unit  area  at  the  specimen-shell  interface. 


The  temperature  drop  across  the  brass  shell  is 

T . . T _ Oaln(c/b) 
rb  rc  v > 

Kg 

where  T is  the  temperature  at  the  outer  surface  of  the  shell,  r * c 


(12) 


and  kg  is  the  thermal  conductivity  of  the  shell  material. 

In  the  apparatus  configuration  shown  in  figure  1,  the  temperatures 
measured  correspond  to  and  Tc  as  defined  above.  Combining  equations 
(6),  and  (9)  through  (12),  we  can  eliminate  all  of  the  temperatures 
involved  except  and  Tc  and  thus  obtain,  instead  of  (6), 


^ Qa  /n(b/a) 

* <Tb.Ta>  = Th  - V;  L1  ’ 


J2sL 


Th  - T- 


.-1 


hSsJlH  + tn(c/b)  + Ra  . Rb  H 

1.  1.  — T —J , . 


(13) 
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The  magnitude  of  the  correction  terms  in  eq.  (13)  is  more  readily  seen 
if  we  define 


app 


Qa  ln(b/a) 
Th  - Tc 


(14) 


as  the  apparent  value  of  thermal  conductivity  obtained  by  neglecting  the 
correction  terms;  combining  (13)  and  (14)  we  obtain 

»*«b.v  - >*.„  [.  - ^ * ¥ * ?)]"  <‘5> 

Equation  (15)  is  the  basic  equation  to  be  used  for  calculating  mean 
thermal  conductivity  values  if  the  heat  is  assumed  to  flow  radially.  The 
quantities  a,  b,  Q,  T^,  and  Tc  involved  in  the  definition  of  \*app 
(eq.  (14))  are  all  directly  measurable.  Ihe  quantities  kc,  ks,  Ra,  and 
Rb  (required  for  the  correction  terms  in  (15))  must  either  be  known  or 
the  corresponding  correction  terms  shown  to  be  negligible. 

4.3.  Longitudinal  Heat  Flow 

The  analysis  given  above  is  based  on  the  assumption  of  purely 
radial  heat  flow  in  the  alumina  core  and  in  the  specimen.  Owing  to  the 
finite  length  of  the  apparatus,  consideration  must  be  given  to  the  effects 
of  longitudinal  heat  flow  in  the  core  and  in  the  specimen. 

Peavy  [2]  has  considered  the  effect  of  longitudinal  heat  flow  in  an 
apparatus  such  as  that  used  in  the  present  investigation.  His  analysis 
is  exact  for  the  case  in  which  both  the  core  and  the  specimen  have  constant 
thermal  conductivity  and  the  heat  generation  per  unit  length  does  not 
vary  along  the  core.  Flynn  [l]  has  considered  the  effect  of  longitudinal 
heat  flow  along  the  core  for  the  case  in  which  the  heat  generation  per 
unit  length  is  a linear  function  of  temperature.  His  analysis  neglects 
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longitudinal  heat  flow  in  the  specimen  and  assumes  a constant  thermal 
conductivity  (i.e.,  independent  of  temperature)  for  the  core. 

In  the  present  investigation  the  very  large  temperature  differences 
within  the  apparatus  make  it  imperative  that  the  temperature  dependence 
of  the  thermal  conductivity  of  the  core,  the  temperature  dependence  of 
the  thermal  conductivity  of  the  specimen,  and  the  temperature  dependence 
of  the  heat  generation  per  unit  length  all  be  taken  into  account.  Peavy's 
analysis  [2]  would  be  very  difficult  to  adapt  to  the  case  in  which  all 
these  parameters  have  large  temperature  dependences.  Flynn’s  analysis  [l] 
can  easily  be  modified  to  include  the  effect  of  the  temperature  dependence 
of  the  thermal  conductivity  of  the  core  and  with  a little  more  effort  be 
extended  to  include  a more  complicated  (than  linear)  temperature  depen- 
dence for  the  heat  generation  per  unit  length  in  the  core.  However,  his 
analysis  is  inherently  incapable  of  dealing  with  the  effects  of  longitudi- 
nal heat  flow  in  the  specimen. 

Since  the  apparatus  design  permits  direct  experimental  measurement 
of  the  temperature  distribution  along  the  axis  of  the  core,  we  have  more 
information  than  is  available  in  boundary  value  problems  in  which  only  the 
conditions  on  the  outside  of  the  system  are  prescribed.  This  additional 
information  enables  us  to  compute  the  effects  of  longitudinal  heat  flow 
in  a rather  straightforward  manner  with  only  minor  approximations  being 
required.  If  necessary  even  these  approximations  could  be  circumvented 
by  a simple  process  of  iteration. 
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Let  us  first  consider  the  situation  which  exists  in  the  alumina  core. 
The  thermal  conductivity  of  high-density,  high-purity  alumina  decreases 
by  a factor  of  more  than  six  between  room  temperature  and  1600  °C.  The 
electrical  resistivity  of  the  607,  Pt  - 407»  Rh  alloy  used  for  the  heater 
winding  increases  by  a factor  of  three  between  room  temperature  and 
1600  °CU  Both  the  thermal  conductivity  of  the  core  and  the  electrical 
resistivity  of  the  heater  (and  hence  the  heat  generation  per  unit  length) 
are  explicit  functions  of  temperature.  Since  we  experimentally  measure 
the  temperature  distribution  along  the  axis  of  the  core  and  since  the 
temperature  at  the  axis  of  the  core  is  essentially  the  same  as  the 
temperature  of  the  heater  winding,  we  can  determine  the  thermal  conductivity 
of  the  core  (inside  the  circle  of  heater  wires)  and  the  heat  generation 
per  unit  length  as  functions  of  position.  For  a given  test  we  can  treat 
both  of  these  quantities  as  if  they  were  explicit  functions  of  position. 

Consider  an  element  of  the  core  of  length  dz.  The  heat  generated 
in  this  element  is  ml23dz,  where  m is  the  number  of  heater  wires,  I is 
the  electric  current  flow  through  the  heater  winding,  and  3 is  the 
electrical  resistance  per  unit  length  of  heater  wire.  The  radial  heat 
flow  from  the  convex  surface  of  this  element  into  the  specimen  is 
2-rraQdz,  where  now  Q = Q(z)  varies  with  longitudinal  position.  The  heat 
flowing  into  this  element  in  the  positive  z-direction  is  -C(dv/dz),  where 
C is  the  ("posi t ion-dependent ")  thermal  conductance  of  the  core, 
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including  beater  wires,  and  v is  the  average  temperature  of  the  core 

at  longitudinal  position  z,  The  beat  leaving  the  element  in  the 

positive  z~direction  is 

-C(dv/dz)  (d/dz) [-C(dv/dz) ] dz . 

Under  steady-state  (i.e„  temperature  independent  of  time)  conditions, 
conservation  of  energy  enables  us  to  write 

(c#)  -2rraQ(z)  = -ml8g(z)  (16) 

dz  \ dz/ 


as  the  differential  equation  relating  T?(z) , Q(z) , and  $(z).  If  v and 
3 are  known  or  measured,  the  heat  flux  into  the  specimen  at  any  longi- 
tudinal position  is 


Mathematically, 
is  given  by 


, s 1 f t2q  m ^ l r dvNi“! 

Q<z>  = 2^  LmI  e + Tz  [c  Tz!  j . 


(17) 


the  average  core  temperature  at  any  longitudinal  position 


v(z)  * “2 
a 


rv(r ,z)dr 


(18) 


With  this  definition  of  v,  eq.  (17)  is  completely  rigorous. 

Except  very  near  the  ends  of  the  apparatus,  most  of  the  heat  flow 
in  the  core  is  radial  rather  than  longitudinal.  Hence,  to  a good  approx- 
imation, the  temperature  at  any  position  in  the  core  is  given  by 

v(0,z)  , 0 < r < r ' 


v(r,z)  =- 


x ml2  8 T 1 dr 

(0,z)  - j , rrr  > r'  < r < 


(19) 


1A 


where  v(0,z)  is  the  experimentally  determined  temperature  distribution 
along  the  axis  of  the  core.  Substituting  this  expression  into  (18) , the 
average  temperature  of  the  core  at  any  cross  section  is  given  very 
closely  by 

V(Z)  = v(0,z)  - rdr  ^ £■  . (20) 

For  the  relatively  small  radial  temperature  differences  between  r = r* 
and  r = a we  can  justifiably  neglect  the  temperature  dependence  of  kc 
and  thereby  reduce  (20)  to 


v (z)  = v(0,z)  - 


tt a~kc  Jr» 


r Xn  r»  dr 


= v(0,z)  - r^n  ^ -X  + ^ 

4rrkc  L r ad 


(21) 


This  equation  can  be  substituted  into  (17)  to  enable  calculation  of 

Q(z)  , the  radial  heat  flux  into  the  specimen  at  all  longitudinal  positions. 

By  writing  a differential  equation  (eq.  (16))  in  terms  of  v,  the 
average  core  temperature,  we  were  able  to  determine  the  heat  flux,  Q(z) , 
into  the  inner  surface,  r = a,  of  the  specimen.  We  now  show  the  analysis 
needed  to  permit  computation  of  the  thermal  conductivity  of  the  specimen. 
Consider  an  apparatus,  having  the  horizontal  cross  section  shown  in 
figure  1,  which  is  of  finite  length,  extending  from  z = 0 to  z = i,  with 
the  flat  ends  water  cooled.  We  define  v = 0 at  r = b,  neglect  contact 
resistances,  and  write  the  boundary  conditions  for  the  specimen  as 
a < r < b z = 0 and  z = i v **  0 i 


r = b 


r = a 


0 < z < H, 

0 < z < i 


v = 0 


(22) 


-x(v)  fr  = Q<z>  • 

or 
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In  order  to  simplify  the  boundary  condition  at  the  surface  of  the  core, 

r = a,  we  introduce  a new  variable  defined  by 

y = 7"*  I X (v)dv  ” \ ^ v (23) 

ko  0 \0  > 

where  \Q  ~ \(v  ~ 0) e In  terms  of  y,  the  boundary  conditions,  (22),  can 
be  written  as 

a < r <;  b z“0  and  z ~ l y = 0 y 

r-b  0 < z <.  $,  y = 0 ? (24) 

r = a 0 < z < y “X0  ^ - Q(z)  • 

With  these  boundary  conditions,  the  solution  of  the  Laplacian  equation  for 
the  potential,  y,  is  given  by  Carslaw  and  Jaeger  [3,  p.  221]  as 


where 


00 

V l K F0(mTr/x;nnb/^)  . nrrz  /,cs 

7it)  * o L mTte;“b/7)  SlnT  » ( ' 


n=l 


-t, 


2 . nrrz  . . 

a„  -7J0'Q<*  > sin~  «*«’ 


(26) 


F0(x;y)  - I0(x)K0(y)  - IQ(y)K0(x)  , (27) 

and  Fx  (x;y)  = Ix  (x)K0(y)  + I0(y)iq  (x)  . (28) 

I„ (x)  and  (x)  being  modified  Bessel  functions  of  the  first  and  second 
kind,  respectively,  and  of  order  n. 

Evaluation  of  (25)  at  r = a and  comparison  with  (23)  enables  us  to 

write 

00 

,*/„>  - - V * k F0(rma/i;nnb/£)  fi.„  nrrz  (29) 

k vv'  v(a,z)  L 


n=l 


nn  F1  (nna/y  ;nnb/j0) 
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Since  we  know  that  the  heat  flow  will  be  essentially  radial  near  the 
center  of  the  apparatus,  it  is  useful  to  rewrite  (29)  in  the  form 


X*(V)  - «?>  «■*"<»»'«>[!  - „(,)]  . 

v(a,z)  J 


(30) 


where  U(z),  a correction  term  for  the  effects  of  longitudinal  heat  flow 
in  the  specimen,  is  given  by 


U(z)  = _L_  A f,  1 l F0(nna/f  ;nnb/Q] 

Q(z)  • V ’ 4n(b/a)  * rma  Fj  (ima/f ,;nnb/£)J 


HIT  7j  / O 1 \ 

in  — . (31) 


In  deriving  eq.  (30),  we  neglected  the  influence  of  contact  resis- 
tances. Since  U(z)  is  a small  correction,  we  assume  that  it  is  still 
valid  even  if  contact  resistances  are  taken  into  account.  Comparing 


eqs.  (14),  (15),  and  (30),  we  write 


* (Tb,Ta)= 


\*app  [1  - U(z>] 


TT 


1 * + JiBlcM  + Ra  + Rb 

N kc  ks  a b 


(32) 


where  \n  now  has  the  more  general  definition 
app 

* Q(z)  ^n(b/a) 


app 


Th  - Tc 


(33) 


where  Q(z)  is  given  by  eq.  (17).  Equation  (32)  is  the  equation  used 
the  present  investigation  to  compute  \ . 


m 
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4.4.  Apparatus 


A vertical  cross  section  of  the  apparatus  designed  and  built  for 
the  thermal  conductance  measurements  is  shown  in  figure  2.  The  specimen 

is  contained  in  the  annular  space  between  the  central  ceramic  core,  which 
is  supported  at  both  ends,  and  the  inner  brass  shell.  Cooling  water 
flows  in  the  3-ram  annulus  between  this  inner  brass  shell  and  a coaxial 
outer  brass  shell. 

The  central  core  is  an  extruded  alumina  rod,  46  cm  long  and  1.25  cm 
in  diameter.  A horizontal  cross-section  of  this  core  is  shown  in  figure  3. 
Sixteen  equally-spaced  holes,  0.09  cm  in  diameter,  extend  the  entire 
length  of  the  rod.  The  centers  of  these  holes  form  a circle  of  0.44  cm 
radius.  The  core  heater,  which  provides  the  heat  flowing  radially  through 
the  specimen,  consists  of  lengths  of  platinum-40  percent  rhodium  (0.06  cm 
diameter)  wire  threaded  back  and  forth  through  the  sixteen  holes.  The 
core  is  held  from  beneath  by  a ceramic  insulating  support  designed  to 
permit  free  expansion  of  the  heater  wire.  The  upper  end  of  the  core 
passes  through  a hole  in  the  removable  flange  at  the  top  of  the  apparatus; 
this  hole  is  a sufficiently  loose  fit  to  permit  free  expansion  of  the 
core0  Current  leads  and  voltage  taps  are  attached  to  the  heater  winding 
at  the  upper  end  of  the  core. 
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Figure  2. 


Vertical  cross  section  of  the  apparatus  for  measuring 
the  thermal  conductance  of  soils. 
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Figure  3. 


Horizontal  cross  section  of  the  ceramic  core. 


In  ,the  center  of  the  alumina  core  there  is  an  axial  hole  0.25  cm  in 
diameter.  This  hole  accommodates  a thermocouple  which  can  be  moved  ver- 
tically by  exterior  manipulation.  This  thermocouple  is  fabricated  from 
0.04  cm  diameter  platinum  and  platinum  - 10  percent  rhodium  wire  and  is 
contained  in  double-bore  alumina  tubing  which  is  a snug  slip  fit  in  the 
thermocouple  well.  The  temperature  measured  at  the  midlength  of  the 
apparatus  is  designated  T^.  This  thermocouple  can  also  be  utilized  to 
obtain  the  longitudinal  temperature  distribution  along  the  core  - this 
information  is  needed  to  make  corrections  for  longitudinal  heat  losses. 

With  the  exception  of  the  ceramic  core,  all  surfaces  in  contact  with 
the  specimen  are  water-cooled  to  maintain  them  at  room  temperature.  As 
shown  in  figure  2,  the  cooling  water  enters  the  center  of  the  bottom  of 
the  apparatus,  passes  upward  in  the  annulus  between  the  inner  and  outer 
brass  shells,  and  exits  at  the  top  of  the  apparatus.  A circulating  system, 
shown  in  figure  4,  is  used  to  maintain  a constant  temperature  in  the 
apparatus.  The  water  leaving  the  apparatus  is  pumped  through  a heat  ex- 
changer which  cools  the  water  to  a temperature  lower  than  that  desired  at 
the  inlet  to  the  apparatus.  From  the  heat  exchanger,  the  water  enters  a 
small  commercial  hot  water  heater  where  it  is  reheated  to  the  desired 
temperature  for  circulation  to  the  apparatus.  A resistance  thermometer 
at  the  exit  from  the  tank  is  used  to  control  the  power  to  the  heater  in 
the  tank  and  thus  maintain  the  desired  water  temperature. 
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Figure  4.  The  circulating  system  used  to  cool  the  apparatus. 
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TWo  copper  versus  constantan  thermocouples  (0.025  cm  diameter)  are 
attached  to  the  outer  surface  of  the  inner  brass  shell  at  the  midplane 
of  the  apparatus.  The  junctions  of  these  thermocouples  are  thermally 
insulated  from  the  cooling  water  by  plastic  electrical  tape.  The  average 
of  these  two  temperatures  is  designated  Tc#  Additional  copper-constantan 
thermocouples  are  used  to  monitor  temperatures  at  different  locations  in 
the  circulating  system. 

Individual  ice  junctions  are  used  with  the  platinum  versus  platinum- 
10  percent  rhodium  thermocouple  used  to  measure  T^.  The  constantan  leads 
of  the  copper  versus  constantan  thermocouples  are  brought  to  an  isothermal 
zone  box  at  room  temperature.  A thermocouple  with  one  junction  in  the 
zone  box  and  one  in  an  ice  bath  is  placed  in  series  with  a double-pole 
selector  switch,  so  that  each  measuring  thermocouple  is  automatically 
referenced  against  the  ice  bath[4]. 

The  ceramic  core  heater,  which  provides  the  heat  flowing  radially 
through  a specimen,  is  fed  a-c  power  from  a saturable  core  reactor  which 
is  controlled  by  a current-adjusting-type  proportional  controller  incorpor- 
ating automatic  rate  and  reset  action.  This  controller  has  been  modified 
to  operate  in  either  of  two  modes.  It  can  adjust  the  power  to  the  core 
heater  so  as  to  maintain  the  temperature,  T^,  at  the  center  of  the  core 
at  a constant  value.  Alternatively,  it  can  control  the  saturable  core 
reactor  so  as  to  provide  constant  current  to  the  core  heater.  The  power 
dissipated  in  the  core  heater  is  determined  by  measuring  the  current  through 
the  heater  and  the  voltage  drop  across  the  heater  utilizing  one-quarter- 
percent-accuracy  electrodynamic  voltmeters  and  ammeters. 


5.  Testing  Procedure  and  Calculation  Procedure 
5*1.  Specimen  Preparation 

Sandia  Corporation  provided  us  with  desired  dry  bulk  density  values 
for  each  of  the  soil  samples  supplied.  These  desired  densities  were 
based  on  the  maximum  densities  attained  in  a series  of  tests  using 
A STM  Method  of  Test  D698.  This  test  method  provides  a specified  compac- 
tion procedure  based  on  dropping  a weight  a fixed  distance  a given  number 
of  times.  This  compaction  procedure  had  been  carried  out  on  each  soil 
for  a number  of  samples  having  different  moisture  contents;  the  moisture 
contents  which  resulted  in  the  highest  dry  bulk  densities  were  reported 
to  us. 

Adaptation  of  ASTM  D698  to  the  geometry  of  our  thermal  conductivity 
apparatus  required  considerable  effort.  We  were  successful  in  developing 
a compaction  procedure,  based  on  ASTM  D698,  which  enabled  rapid,  reliable 
compaction  of  the  soil  samples  with  reproducible  densities.  The  bulk 
densities  obtained  with  our  compaction  rig  did  not  always  agree  with 
those  densities  obtained  with  ASTM  D698  for  similar  samples  at  nominally 
the  same  moisture  content  and  compaction  energy.  We  were  always  able 
to  come  reasonably  close  to  the  desired  bulk  densities  with  only  a few 
trials,  however. 
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Each  thermal  conductivity  sample  was  compacted  into  a new  brass  can 
(the  inner  brass  shell  shown  in  figure  2).  This  can  had  a centered 
hole  in  the  bottom  which  held  the  lower  end  of  the  ceramic  core  (heater) 
coaxially  within  the  can.  During  the  compaction  operation , a combination 
weight  guide  and  core  support,  as  shown  in  figure  5,  was  affixed  to  the 
upper  end  of  the  brass  can  which  was  to  contain  the  specimen.  The  core 
support  held  the  upper  end  of  the  ceramic  core  (heater)  coaxially  within 
the  can  during  compaction;  the  core  support  was  spring-loaded  to  hold 
the  ceramic  core  down  during  compaction.  A brass  weight  was  used  to 
compact  the  soil  in  the  specimen  container.  This  weight  had  an  outside 
diameter  slightly  smaller  than  the  inside  diameter  of  the  specimen 
container  and  was  provided  with  a keyhole-shaped  slot  to  clear  the  core 
support  and  web  shown  in  figure  5.  An  attachment  extending  above  the 
top  of  the  weight  prevented  the  weight  from  rotating  even  when  the 
weight  was  totally  within  the  specimen  can.  The  weight  guide  could  be 
rotated  when  desired  so  that  the  keyhole-shaped  slot  in  the  weight 
would  not  always  impinge  on  the  same  position  of  the  specimen. 

Specimen  preparation  consisted  of  the  following  steps: 

1.  A small  amount  of  the  soil  was  oven-dried  to  constant  weight 
at  110  °C  in  order  to  determine  the  existing  moisture  content. 
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re  5.  Rig  used  for  compacting  the  soil  samples. 
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2. 


Meanwhile  the  brass  can  (which  would  later  receive  the  specimen) 
and  the  ceramic  core  (heater)  were  weighed. 

Sufficient  additional  moisture  was  added  to  a soil  sample  to 
attain  the  moisture  content  desired. 

4.  The  compaction  rig  shown  in  figure  5 was  affixed  to  the  top 
of  the  specimen  container  with  the  heater  centered  in  place. 

5.  A mass  of  moisture-conditioned  soil  equal  to  one-tenth  of  the  total 
mass  needed  to  achieve  the  desired  density  was  placed  in  the 
specimen  container. 

6.  The  weight  was  dropped  a fixed  distance  a given  number  of  times 
to  achieve  the  desired  compaction  energy. 

7.  Another  layer  of  soil  was  placed  in  the  specimen  container  and 
this  layer  was  compacted  as  above.  This  procedure  was  repeated 
until  the  specimen  container  was  filled. 

8.  The  filled  specimen  container  was  oven-dried  to  constant  weight 
at  110  °C. 

9.  The  specimen  and  container  were  weighed  and  the  specimen  dry 
bulk  density  was  computed. 

10.  If  the  dry  bulk  density  obtained  was  sufficiently  close  to  the 
desired  value,  testing  proceeded.  Otherwise,  the  specimen  was 
removed  and  re-packed  using  a different  compaction  energy. 
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A few  of  the  soil  samples  did  not  lend  themselves  to  the  above 
compaction  procedure . The  specimens  of  Ottawa  sand,  dune  sand,  and 
coastal  plains  clay  were  packed  to  the  desired  density  by  filling  the 
specimen  container  slowly  with  dry  soil  while  rapping  it  with  a small 
hammer. 

5.2.  Testing  Procedure 

A specimen  container,  with  the  ceramic  core  heater  and  the  canpacted 
specimen  installed,  was  placed  in  the  outer  brass  shell  as  shown  in 
figure  2.  Rubber  O-rings,  not  shown  in  figure  2,  provided  water  tight 
seals  at  the  top  of  the  apparatus.  The  current  leads  and  voltage  taps 
were  attached  to  the  heater  leads  extending  from  the  top  of  the  ceramic 
core.  A thermocouple  was  installed  in  the  axial  hole  in  the  ceramic 
core  so  that  the  hot  junction  was  at  the  mid-height  of  the  apparatus. 

The  upper  end  of  this  thermocouple  was  attached  to  a height  gage 
which  permitted  moving  the  thermocouple  vertically  by  known  distances. 

The  circulating  system  shown  in  figure  4 was  activated  to  maintain 
the  outside  of  the  specimen  container  at  the  desired  temperature.  This 
involved  adjusting  the  amount  of  chilled  water  passing  through  the  heat 
exchanger  as  well  as  the  setting  of  the  temperature  controller  connected 
to  the  resistance  bulb  in  the  water  tank. 
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The  specimen  heater  was  energized  with  the  proportional  controller 
operating  in  a cons tant- temperature  mode  so  as  to  bring  and  maintain 
the  temperature  at  the  core  thermocouple  to  the  desired  value.  This 
procedure  effectively  prevented  significant  temperature  overshoot, 
allowing  monotonic  heating  of  the  hot  face  of  the  sample  but  yet  per- 
mitting rapid  equilibrium  to  be  attained.  When  the  temperature  indicated 
by  the  core  thermocouple  was  steady  and  when  the  electrical  power 
required  to  maintain  this  temperature  had  also  leveled  out  to  a steady 
value,  the  proportional  controller  was  switched  to  a cons tant -current 
mode  so  as  to  maintain  the  current  to  the  heater  at  the  value  it  had  had 
before  switching  over. 

After  any  small  disturbances  due  to  mode-changing  had  had  time  to 
die  out,  several  successive  readings  were  taken  of  the  current  through 
the  heater,  the  voltage  drop  across  the  heater,  the  emf  of  the  Pt  - 107c  Rh 
Pt  thermocouple  (at  the  mid-height  of  the  apparatus)  in  the  axial  hole  in 
the  ceramic  core,  and  the  emfs  of  the  two  copper jconstantan  thermo- 
couples on  the  outside  of  the  specimen  container.  After  satisfactory 
data  were  obtained  with  the  core  thermocouple  at  the  mid-height  of  the 
apparatus,  the  thermocouple  junction  was  moved  to  a position  approximately 
level  with  the  bottom  of  the  specimen.  The  emf  of  the  core  thermocouple 
was  read  at  this  position  and  then  at  approximately  2.5  cm  intervals 
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along  the  entire  length  of  the  specimen,  time  being  allowed  at  each 
position  for  the  thermocouple  junction  to  reach  thermal  equilibrium  with 
the  ceramic  core.  The  core  thermocouple  junction  was  then  returned  to 
the  mid-height  of  the  apparatus  and  readings  were  again  taken  of  the  heater 
current,  heater  voltage,  and  thermocouple  emfs. 

Following  completion  of  one  set  of  data,  the  heater  controller  was 
switched  to  constant  temperature  mode  and  the  heater  and  both  controllers 
were  adjusted  to  bring  the  specimen  to  the  next  condition  desired.  Data 
were  normally  taken  at  200  °C  intervals  to  about  1100  °C  (hot-side)  and 
then  at  100  °C  intervals  until  either  the  heater  burned  out  or  the  core 
thermocouple  ceased  to  function  properly. 

For  the  first  tests  conducted,  the  heater  was  fabricated  by  connecting 
all  sixteen  heater  wires  in  series.  With  this  wiring  arrangement,  the 
current  leads  to  and  from  the  heater  were  in  adjacent  holes  in  the 
ceramic  core.  This  resulted  in  large  voltage  gradients  within  the 
alumina  core  with  consequent  dielectric  breakdown  of  the  alumina  at 
core  temperatures  around  1400  °C.  The  heater  wiring  was  subsequently 
changed  to  a series-parallel  arrangement  in  which  the  heater  was  wired 
as  two  sets  of  eight  series-connected  wires  with  these  two  sets  con- 
nected in  parallel.  This  greatly  reduced  the  voltage  gradients  in  the 
ceramic  core  and  we  had  no  further  trouble  with  dielectric 
breakdown. 


30 


Another  problem  that  caused  some  concern  was  that  of  a-c  pickup 
at  high  temperatures  by  the  core  thermocouple.  The  physical  arrangement 
of  the  apparatus  made  it  difficult  to  shield  this  thermocouple.  Pro- 
bably one  of  the  best  procedures  would  have  been  to  use  commercially- 
available  me tal-sheathed  thermocouples  fabricated  by  swaging  techniques. 

The  temperatures  involved  would  have  required  the  use  of  thermocouples 
with  noble-metal  sheaths,  thus  greatly  increasing  the  cost  of  the 
thermocouples.  Since  it  was  desired  to  use  a new  thermocouple  for  each 
specimen  (the  high  temperatures  could  cause  thermocouples  to  lose  their 
calibration),  the  total  cost  would  have  been  quite  high.  It  was  found 
possible  to  use  unshielded  thermocouples  by  reading  the  output  using  a 
electronic  integrating  digital  voltmeter  having  140  dB  common-mode 
rejection. 

5.3*  Calculation  Procedure 

A computer  program  was  written  to  enable  calculation  of  thermal 
conductivity  values  from  the  data  which  were  obtained.  The  input  to 
the  computer  consisted  of  the  voltmeter  and  ammeter  readings  and  the 
thermocouple  positions  and  emfs  for  each  test.  The  calculation  procedure 
was  as  follows: 

1.  The  ammeter  and  voltmeter  readings  were  adjusted  slightly,  if 

necessary,  by  sub-programs  which  contained  the  meter  calibrations . 
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2.  Temperatures  were  computed  from  the  thermocouple  emfs  using 
sub-programs  which  contained  the  thermocouple  calibrations. 

3«  The  power  output  per  unit  length  in  the  absence  of  longitudinal 
heat  losses  was  computed  for  each  position  along  the  ceramic 
core  using  the  leading  term  in  eq.  (17).  The  resistance  per 
unit  length,  3,  of  the  heater  wire  was  computed  using  an 
equation  representing  NBS-obtained  data  on  the  temperature- 
dependence  of  Pt  - 40%  Rh  alloy.  It  was  assumed  that  the 
heater  wire  was  at  the  temperature,  T^,  measured  at  the 
center  of  the  core. 

4.  The  average  core  temperature  at  each  longitudinal  position  was 
computed  using  eq.  (21). 

5.  A fourth-order  polynomial  was  fitted  (by  least-squares)  to  the 
(usually  19)  values  of  v(z).  From  this  equation  and  an  equation 
which  was  used  to  compute  C,  the  longitudinal  thermal  conduc- 
tance of  the  core,  the  second  term  in  eq.  (17)  was  evaluated. 

6.  From  these  values  of  Q(z)  , the  mean  thermal  conductivity, 

X (Tfc,Ta),  was  computed  at  each  position  using  eqs.  (32)  and 
(33)*  The  computer  program  contained  a value  for  kg , the 
thermal  conductivity  of  the  brass  shell,  and  an  equation  for 
computing  , the  thermal  conductivity  of  the  ceramic  core,  as 
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a function  of  temperature,,  Since  no  knowledge  was  available 
regarding  the  thermal  contact  resistances,  Rfl  and  R^  were 
set  equal  to  zero.  A separate  subprogram  was  used  to  compute 
U(z),  the  correction  for  longitudinal  heat  flow  in  the  specimen. 

7.  An  equation  of  the  form  \ (25, T ) = a + bT  + cT  * + dT  n (n  3s  3) 

a a a a 

* 

was  fitted  (least-squares)  to  the  computed  X (T.  ,T  ) values 

d a 

corresponding  to  the  mid-height  of  the  apparatus,  the  X*  ♦•lues 
being  adjusted  slightly  so  that  they  all  corresponded  to  a 
cold  side  temperature  of  = 25  °C. 

8.  The  thermal  conductivity  was  computed  from  this  equation  using 
the  relation 

\(Ta)  = \*(T5.Ta>  + (Ta  - Il)5X  (TbtV  (34) 

dTa 

which  is  obtained  by  differentiation  of  eq.  (5). 

it  it 

9.  The  computed  \ values  and  the  fitted  curves  for  X and  \ were 
stored  on  a magnetic  tape  which  later  was  used  to  generate  the 
plots  (figures  6 through  16)  in  the  next  section. 
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6.  Results 


Prior  to  undertaking  measurements  of  the  thermal  conductivity  of 

any  of  the  supplied  soil  samples , we  measured  the  thermal  conductivity 

of  a sample  of  diatomaceous  earth.  We  had  previously  [l]  measured  the 
thermal  conductivity  of  this  same  lot  of  diatomaceous  earth  in  an  appa- 
ratus in  which  both  the  hot  and  the  cold  sides  of  the  specimen  were 
heated  so  that  X could  be  measured  directly  rather  than  being  computed 
from  X* . Further,  in  this  earlier  apparatus  the  ends  were  heated  so 
that  there  was  no  need  to  make  corrections  for  longitudinal  heat  flow 
in  the  core  or  in  the  specimen.  Measurement  of  this  lot  of  diatomaceous 
earth  in  the  present  apparatus  provided  a check  on  the  corrections  for 
longitudinal  heat  flow  (which  were  rather  large  because  of  the  low 
thermal  conductivity  of  diatomaceous  earth)  and  also  provided  a check 
on  the  procedure  used  to  compute  X from  the  measured  \ values. 

JL 

The  solid  triangles  shown  in  figure  6 represent  the  \ (25, T«)  data 
points  obtained  in  the  present  apparatus  plotted  versus  Ta,  the  hot 
side  temperature.  The  solid  curve  is  the  cubic  of  least-squares  fit 
to  these  data.  The  dashed  curve  represents  the  X(Ta) -values  computed 
using  eq.  (34).  The  open  circles  represent  the  smoothed  thermal  con- 
ductivity values  obtained  earlier  by  Flynn[l]  on  diatomaceous  earth  from 
the  same  lot  as  was  used  in  the  present  measurements.  The  agreement 
is  seen  to  be  good. 
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Figures  7 through  16  show  the  measured  X*-values,  the  least-squares 
X -curves,  and  the  derived  ^-curves  for  each  of  the  nine  soil  samples 
supplied  by  Sandia0  The  data  correspond  to  monotonic  heating  to  the 
maximum  temperatures  shown.  Measurements  have  been  completed  on  two 
samples,  designated  by  solid  triangles  and  open  squares,  at  nominally  the 
same  density,  for  calcareous  soil,  dune  sane,  and  coastal  plains  clay. 

At  the  request  of  Sandia,  the  second  run  on  Ottawa  sand  was  made  with  a 
lower  density  than  was  used  for  the  first  specimen. 

Figure  17  shows  \ for  all  of  the  specimens  measured  to  date  while 
figures  18  and  19  illustrate  the  corresponding  values  of  Note  that 
each  set  of  curves  tends  to  converge  to  a common  region  at  high  tempera- 
tures. The  key  to  these  curves  is  given  below. 


Symbol 

Specimen 

Average  Density 

kg  m“3 

DE 

diatomaceous  earth 

150 

C 

calcareous 

2000 

GD 

granitic  detrital 

1920 

DS 

dune  sand 

1570 

M 

magnesian 

1760 

P 

podzol 

1760 

CP 

coastal  plains 

1335 

L 

laterite 

1490 

EP 

estancia  playa 

1540 

OS, 

Ottawa  sand 

1760 

OSn 

Ottawa  sand 

1570 
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TEMPERATURE  . °C 


Figure  6.  Thermal  conductivity  of  diatomaceous  earth  at  a density  of 

0.15  x 103kg  m"3.  The  solid  curve  represents  an  equation  of 

the  form  X*(25,T  ) = a + bT  + cT  2 + dT  3.  The  open  circles 
a a a a 

represent  the  smoothed  values  obtained  earlier  by  Flynn[l]. 


THERMAL  CONDUCTIVITY. 


I 


Figure  7.  Thermal  conductivity  of  calcareous  soil  at  densities  of  2.01  x 10° 

and  1.99  x 103  kg  m~3.  The* solid  curve  represents  an  equation 

of  the  form  \*(25,T  ) = a + bT  + cT  2 + dT  13. 

v a'  a a a 


37 


THERMAL  CONDUCTIVITY 


Figure  8.  Thermal  conductivity  of  granitic  detrital  soil  at  a density 


of  1*92  x 103  kg 


of  the  form  A*(25,T  ) = 

a 


m w.  The  solid  curve  represents  an  equation 
a + bT  + cT  2 + dT  3. 
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THERMAL  CONDUCTIVITY. 
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Figure  9.  Thermal  conductivity  of  two  specimens  of  dune  sand,  each  at  a 
density  of  1.57  x 103  kg  m"3 . The  solid  curve  represents  an 
equation  of  the  form  X*(25,T  ) = a + bT  + cT  dT  6. 

3 Si  3 3 
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c n o j i u u y 3 h i 


Figure  10.  Thermal  conductivity  of  magnesian  soil  at  a density  of 

1.76  x 10s  kg  m~3.  The  solid  curve  represents  an  equation 

of  the  form  \*(25,T  ) = a + bT  + cT  2 + dT  3. 

v a'  a a a 
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Figure  11. 


Thermal  conductivity  of  podzol  soil  at  a density  of  1.76  x 10s 
kg  m“3.  The  solid  curve  represents  an  equation  of  the  form 
X*(25,Ta)  * a + bTa  + cTa2  + dTa4 . 
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THERMAL  CONDUCTIVITY, 
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Figure  12,  Thermal  conductivity  of  coastal  plains  clay  at  densities  of 

1.32  x 103  and  1,35  x 103  kg  m~3.  The  solid  curve  represents 

an  equation  of  the  form  \*(25,T  ) = a + bT  + cT  2 + dT  12. 

a'  a a a 
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Figure  13  • Thermal  conductivity  of  laterite  soil  at  a density  of 

1.49  x 103  kg  m“3.  The  solid  curve  represents  an  equation 

of  the  form  \*(25,T  ) = a + bT  + cT  2 + dT  6 . 

a'  a a a 
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Figure  14.  Thermal  conductivity  of  estancia  playa  soil  at  a density  of 

1.54  x 103  kg t m The  solid  curve  represents  an  equation 

of  the  form  \,f(25,T  ) - a + bT  + cT  2 + dT  11  . 

a a a a 
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Figure  15.  Thermal  conductivity  of  Ottawa  sand  at  a bulk  density  of 

1.76  x 103  kg  m"3.  The  solid  curve  represents  an  equation 

of  the  form  \'V(25,T  ) = a + bT  + cT  2 + dT  3 . 

a a a a 
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THERMAL  CONDUCTIVITY, 


Figure  16.  Thermal  conductivity  of  Ottawa  sand  at  a bulk  density  of 

1.57  x 103  kg  m’3.  The  solid  curve  represents  an  equation 

of  the  form  \*(25,T  ) = a + bT  + cT  2 + dT  3. 

N * a'  a a a 
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Figure  17.  The  effective  thermal  conductivity,  \*(25,Ta),  of  all  of 
the  specimens  measured  to  date  plotted  versus  the  hot- 
side  temperatures. 
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Figure 


18.  The  derived  thermal  conductivity,  \(Ta),  of  all  of  the 

specimens  measured  to  date  plotted  versus  temperature.  The 
curves  are  truncated  near  \ - 1.4  W m”1  °C“1  . 
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F igure 


19.  The  derived  thermal  conductivity,  \(Ta),  of  all  the  specimens 
measured  to  date  plotted  versus  temperature.  The  vertical 
scale  has  been  changed  to  enable  showing  the  high-temperature 
behavior . 
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Appendix  A:  Soil  Sample  Characterization 


A private  testing  laboratory  carried  out  measurements  of  the  parti- 
cle size  distribution  of  all  of  the  soil  samples  supplied  to  NBS.  We 
hope  to  include  these  measurements  in  the  final  report. 

The  University  of  New  Mexico  carried  out  X-ray  diffraction  analyses 
of  the  soil  samples  and  then  made  estimates  of  the  composition  and  melting 
range  of  each  soilfjJ.  These  results  are  presented  below: 


CALCAREOUS  Common  to  arid  and  semi-arid  regions  of  world.  CaO,  one 

of  12  most  plentiful  oxides  of  the  earth’s  crust,  tend?  to 
have  very  high  melting  range  and  will  cause  gassing. 


Soil  Constituents 

Estimated  Amount 

7o 

Calcite 

Montmoril Ionite 

Quartz 

Dolomite 

Muscovite- il lite 
Goethi te 
01 igoclase 
Trace  Minerals 

75 

8 

8 

3 

2 

1 

1 

2 

100 

Estimated  Soil-Melt  Range:  1950  to  2050  C 
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GRANITIC  PETRI TAI  These  soils  with  a high-alkali  feldspar  content  have 

a low-melting  range.  They  are  typical  of  many 
detrital  soils,  which  comprise  about  a quarter  of 
total  land  area  of  world. 


Soil  Constituents 

Estimated  Amount 

% 

Quartz 

45 

Albite 

20 

Microcl ine 

10 

Calcite 

6 

Montmoril Ionite 

7 

Illite-mica 

3 

Chlorite 

3 

Mixed-layer  clays 

4 

Trace  Minerals 

2 

100 

Estimated  Soil-Melt  Range:  1200  to  1320  °C 


DUNE  SAND  Typical  of  barchan  sand  dunes  of  world,  representative  of 
about  77o  of  total  land  area  of  world.  Consists  largely  of 
siliceous  sand,  Si02 , comprising  nearly  60%  of  oxides  of 
earth's  crust. 


Soil  Constituents 

Estimated  Amount 

% 

Quartz 

60 

Oligoclase 

25 

Microcline 

10 

Clinopyroxene 

5 

100 

Estimated  Soil-Melt  Range:  1300  to  1430  °C 
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MAGNESIAN  Representative  of  soils  derived  from  gabbros,  basalts,  and 
dolomites,,  MgO  is  one  of  the  12  most  plentiful  oxides  of 
the  earth’s  crust. 


Soil  Constituents 

Estimated  Amount 

<V 

/o 

Talc 

35 

Penn ini te 

30 

01 igoclase 

25 

CaLcite 

5 

Trace  minerals 

5 

100 

Estimated  Soil-Melt  Range:  1350  to  1435  °C 


PODZOL  Leached  organic  soils  of  woodland  regions  of  temperate  zones 
of  world,  which  comprise  about  a quarter  of  total  land  area 
of  world. 


Soil  Constituents 

Estimated  Amount 

% 

Quartz 

60 

Microcline 

15 

Albi te 

10 

Chlorite 

5 

Hornblende 

5 

Trace  Minerals 

5 

100 

— . — — j 

Estimated  Soil-Melt  Range:  1350  to  1450  cC 
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COASTAL 


LATE RITE 


PLAINS  CLAY  Typical  of  well-weathered  soils,  representative  of 
about  15%  of  total  land  area  of  world.  Soil  may 
seal  and  build  gas  pressure  when  saturated. 


Soil  Constituents 

Estimated  Amount 

7 

/o 

45 

25 

25 

5 

100 


Estimated  Soil-Melt  Range:  1600  to  1620  °C 


Quartz 
Kaolinite 
Montmoril Ionite 
Trace  minerals 


Representative  of  tropical  humid  regions  of  world,  comprising 
more  than  10%  of  total  land  area  of  world.  FeO  and  Fe203 
are  two  of  most  plentiful  oxides  of  the  earth's  crust. 


Soil  Constituents 

Estimated  Amount 

% 

Goethite 

55 

Gibbs ite 

30 

Magnet j te 

5 

Maghemite 

5 

Ilmeni te 

5 

100 

Estimated  Soil-Melt  Range:  1900  to  1950  °C 
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ESTANCIA  FLAYA  Typical  of  very  saline  bolson  playas  of  the  world;  most 

extreme  chemical  corrosive  action  upon  heating* 


Soil  Constituents 

Estimated  Amount 

% 

Gypsum 

31.5 

Calcite 

13.0 

Dolomite 

24.5 

Quartz 

7.5 

Feldspar 

5.0 

Halite 

7.5 

Clays 

5.0 

Trace  minerals 

6.0 

100.0 

Estimated  Soil-Melt  Range:  1350  to  1400  °C 


OTTAWA  SAND  Commercially  available  material  that  might  possibly  be  used 

as  a standard  in  comparing  results  from  different  laboratories. 

Essentially  pure  quartz 


Estimated  Melting  Point:  1720  °C 
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8.  Appendix  B:  Results  of  Differential  Thermal  Analyses  of  Soil  Samples 


The  soil  samples  were  ground  to  pass  through  a 100-mesh  sieve0 
After  drying  each  sample  at  140  °C  for  two  hours,  the  analyses  were  per- 
formed using  an  R.  L.  Stone  #230  differential  thermal  analysis  (DTA) 
instrument.  Each  sample  was  placed  in  a platinum  cup  and  tapped 
lightly  to  fill  possible  air  spaces.  The  aluminum  oxide  powder  used  as 
a reference  material  was  placed  in  a similar  cup.  The  sample  and  refer- 
ence cups  were  installed  within  a platinum  enclosure  in  the  furnace 
chamber.  The  temperature  difference  between  the  sample  and  reference 
cups  was  monitored  by  a differential  thermocouple  with  a junction  in 
each  cup.  The  enclosure  temperature  was  monitored  by  a thermocouple 
having  its  hot  junction  located  between  the  sample  and  reference  cups. 

For  each  soil  sample,  DTA  curves  were  obtained  at  heating  rates  of 
10,  20,  and  30  °C/minute  with  helium  flowing  through  the  sample  enclo- 
sure; on  certain  samples,  cooling  curves  at  one  or  more  rates  were  also 
obtained.  Heating  curves  at  30  °C/rainute  with  dry  air  flowing  through 
the  sample  enclosure  were  also  obtained  for  each  sample.  No  significant 
differences  were  apparent  between  the  DTA  curves  in  helium  and  in  air. 

The  following  table  summarizes  the  DTA  results  obtained  on  heating. 
The  characteristic  temperatures  A,  B,  and  D were  determined  as  shown  in 
figure  20o 
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In  addition  to  carrying  out  differential  thermal  analyses  (DTA)  on 
each  sample  to  about  1400  °C,  thermogravimetric  analyses  (TGA)  were  made 
on  each  sample  to  500  °C  or  higher  in  order  to  determine  the  amount  of 
free  or  bound  water  coming  out  of  the  sample.  The  TGA  measurements 
were  made  on  undried  samples  in  air  at  a heating  rate  of  15  °C/minute. 

Representative  results  of  the  DTA  and  TGA  tests  are  shown  in  figures 
21  through  38 c 


D 


Figure  20.  Scheme  utilized  to  determine  the  temperatures  A,  B,  and  D 

which  characterize  a particular  endotherm  or  exotherm.  Tie 
ordinate  represents  the  temperature  difference  between  the 
sample  and  reference  cups;  the  abscissa  represents  the  tem- 
perature of  the  sample,  enclosure. 
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Differential  Thermal  Analysis  Results  from  Heating  Curves 


Material 

Run 

Atmosphere 

Heating 

Rate 

Endotherm 

Exotherm 

A 

B 

D 

D 

°C/min. 

°C 

°C 

°C 

°C 

Calcareous 

1 

Helium 

20 

639 

729 

826 

1268 

2 

Helium 

10 

— 

— 

(372) 

.... 

3 

Helium 

10 

704 

715 

782 

— 

4 

Helium 

30 

656 

683 

784 

1265 

5 

Helium 

30 

— 

... 

886 

— 

6 

Air 

30 

683 

774 

891 

1262 

Granitic 

1 

Helium 

20 

563 

563 

564 

892 

Detr ital 

683 

683 

738 

— 

— 

— 

1156 

... 

2 

Helium 

10 

560 

560 

563 

... 

— 

— 

720 

— 

3 

Helium 

30 

561 

561 

568 

904 

725 

725 

782 

— 

1119 

1119 

1160 

— 

4 

Air 

30 

562 

562 

570 
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Figure  21.  Thermogravimetric  analysis  of  calcareous  soil. 
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Figure  22.  Differential  thermal  analyses  of  calcareous  soil. 
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Figure  23  Thermogravimetrie  analysis  of  granitic  detrital  soil. 
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Figure  24  Differential  thermal  analyses  of  granitic  detrital  soil. 
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Figure  25  Thermogravimetric  analysis  of  dune  sand. 
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Figure  26  . Differential  thermal  analyses  of  dune  sand. 
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Figure  27.  Thermogravimetric  analysis  of  magnesian  soil. 
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Figure  28.  Differentil  thermal  analysesof  magnesian  soil. 


67 


% ‘39NVH9  1H9I3M 


o 

O 


Figure  29.  Thermogravimetric  analysis  of  podzol  soil. 
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Figure  30.  Differential  thermal  analyses  of  podzol  soil. 
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Figure  31.  Thermogravimetric  analysis  of  coastal  plains  clay. 
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Figure  32.  Differential  thermal  analyses  of  coastal  plains  clay. 
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Figure  33.  Thermogravimetric  analysis  of  laterite  soil. 
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Figure  34.  Differential  thermal  analyses  of  laterite  soil. 
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Figure  35.  Thermogravimetric  analysis  of  estancia  playa  soil  . 


74 


ESTANCIA  PLAYA  SOIL 


i °c  i °c 


0 200  400  600  800  1000  1200  1400 

TEMPERATURE,  °C 


Figure  36.  Differential  thermal  analyses  of  estancia  playa  soil. 
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Figure  37.  Thermogravimetric  analysis  of  Ottawa  sand. 
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Figure  38.  Differential  thermal  analyses  of  Ottawa  sand. 
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